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Decay Rates of Various Bottomonium systems* 

Seyong Kim^ 

High Energy Physics Division, Argonne National Laboratory, 9700 South Cass Avenue, 
Argonne, IL 60439, USA 

Using the Bodwin-Braaten-Lepage factorization theorem in heavy quarkonium decay and production processes, 
we calculated matrix elements associated with S- and P-wave bottomonium decays via lattice QCD simulation 
methods. In this work, we report preliminary results on the operator matching between the lattice expression and 
the continuum expression at one loop level. Phenomenological implications are discussed using these preliminary 
MS matrix elements. 
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Heavy quarkonium decay rates can be written 
as a sum of factorized products of perturbative 
coefficients and non-perturbative matrix elements 
[1], 

T(H -> LH) = 



for the decays into light hadrons and 

T(H -> EM) = 
2Im/,f M (A) 



M a 



H\O* M (A)\0){0\O' n EM (A)\H), 

(2) 



for the electromagnetic decays, where A is a fac- 
torization scale and /„ 's are the coefficients from 
four-fermi vertices. 

For the case of bottomonium decays, we have 
calculated the lattice regulated matrix elements 
in this factorization theorem [2], using the leading 
order non-relativistic QCD lagrangian 



f)2 f)2 

C = 4>UD t U + xHD t + W. 

v 2M' v 2M' 



(3) 



The matrix elements in which we are interested 
are : 

Gi = ( 1 S\^xx^\ 1 S)/M 2 (4) 
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F 1 = (i5|^ x | )(0|x t (-^ D) 2 ^\ 1 S)/M i (5) 

E x = CP\^(i/2) D XX t («72) D ^\ 1 P)/M 4 (6) 

H 8 = { 1 P\^T a xx^T a iP\ 1 P)/M 2 . (7) 

To make use of lattice regulated matrix ele- 
ments (which we have calculated already in [2]) 
for phenomenology, we need to perform operator 
"renormalization" which translates lattice regu- 
lated matrix elements into continuum MS matrix 
elements, 



Of 



MS 
j ' 



(8) 



where Cs are ordered following the velocity scal- 
ing law [3] . Once we calculate the matrix Z, \j , we 

obtain Of^ s by inverting the Z matrix : 



(9) 



From G\, F\, H\, and Hg, we note that neces- 
sary Cs are : 



I = ^xxH, 

Z? 2 = ^ x |0)( | x t(_lS)2^ 

for the S-wave case and 

h x = ^(i/2) DxxH*/2) Dip, 

h 8 = ^T a xx^T a i>, 



(10) 
(11) 

(12) 
(13) 



for P-wave. Here, we dicuss our preliminary res- 
ults on the renormalization of these operators (full 
details will be reported elsewhere [4]). 
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Firstly, we derive Feynman rules for our lat- 
tice lagrangian upto 0(g 2 ) from the lattice quark 
propagator, 

G(x,t+1) = (1 - H /2n) n U\x,t){l - H /2n) 

+ ^,0^ + 1,0' 

with G{x,t) = for t < and H = -A 2 /2M - 
E su b (A is lattice covariant derivative, E su b = 
3(1 - «o)/M , mo = (0\±TrU plaq \0)i,n = 2). 
Then, using these Feynman rules, we calculate 
one loop renormalization of the operators, (10) - 
(13). 

Necessary Feynman diagrams for the part of 
D 2 which contributes to the mixing with the con- 
tinuum operator / are given in Fig 1. A dotted 
line means the temporal component of the gluon 
field and a wavy line means the spatial component 
of gluon field (in the following figures, we use the 
same convention), a) has the lattice V 2 vertex, 
b) has the g 2 A 2 vertex, and c) is the tadpole con- 
tribution. After finding analytic expressions for 
these diagrams, we take p (external momentum) 
— > (except the tadpole contribution) to pick out 
pieces for the continuum I. 




Figure 1 . Feynman diagram for the power diver- 
gent part of D 2 renormalization. 



Similarly, Fig 2 is the Feynman diagrams which 
need to be evaluated for the finite renormalization 
of D 2 operator, a) and b) have the lattice V 2 ver- 
tex, c) has x^i 1 vertex, d) is the tadpole contri- 
bution, e) has the gA ■ V vertex, and f) has the 
g 2 A 2 vertex. We take (analytic expression 

for these diagrams) |^ =0 in order to get the finite 
renormalization of D 2 except for the case of the 
tadpole contribution. 




Figure 2. Feynman diagram for the multiplicative 
renormalization of D 2 operator. 



The finite renormalization of the operator, I, is 
given by the Feynman diagrams of Fig 3. Both a) 
and b) have x^i 1 vertex. Again, to get the finite 
renormalization of I, we set p = in the analytic 
expression for the Feynman diagrams. Expansion 
of these Feynman digrams in terms of the external 
momentum tells us that this operator mixes with 
continuum V 2 at order 0(a s v 2 ) (actually, both 
lattice D 2 and / mix with all continuum terms in 
arbitrary powers oft; 2 ). However, since a s ~ v 2 , 
this effect is in higher order of a s and we neglect 
such mixing effect to be consistent with the order 
we are calculating. 




Figure 3. Feynman diagram for the multiplicative 
renormalization of / operator. 



Fig 4 shows the Feynman diagrams for the fi- 
nite renormalization of the operator, hg. The ver- 
tex factor is T a (g) T a . Expansion in terms of 
the external momentum shows that this term does 
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Figure 4. Feynman diagrams the multiplicative 
renormalization of hg 




Table 1 



numerical values of Feynman integrals 



diagram 


Vegas output 


fig la 


0.4781(8) 


fig lb 


0.7008(2) 


fig lc 


-0.8448(100) 


fig 2a 


0.12476(3) 


fig 2b 


-0.18774(3) 


fig 2c 


0.005360(4) 


fig 2d 


0.1411(100) 


fig 2e 


-0.27241(5) 


fig 2f 


-0.116804(7) 


fig 3 


0.010717(2) 


fig 4a+4b 


-0.001346(3) 


fig 5a+5b 


-0.11773(3) 


fig 6a+6b+6c+6d 


-0.0060293(9) 


(fig 6a+6b+6c+6d) log 


0.0013 



Figure 5. Feynman diagrams for the multiplicat- 
ive renormalization of hi 



not mix with hi. With p = 0, we obtain the fi- 
nite renormalization contribution from these dia- 
grams. 

With p = p 1 = 0, the Feynman diagrams in Fig 
5 gives the contribution to the finite renormaliz- 
ation of the operator, hi. The vertex factor is 
lattice version of V ® V. 




Figure 6. Feynman diagrams for the mixing con- 
tribution from hi to h% 



Finally, Fig 6 shows the Feynman diagrams for 
the mixing of the operator, hi, with hg. The ver- 
tex factor is the same as Fig 5. With zero external 



momentum, these diagrams have a logarithmic in- 
frared divergence. Therefore, these digrams need 
extra care. We subtract logarithmic divergent 
piece from the analytic expression of these dia- 
grams if JV p'j < ir 2 when the integral is evalu- 
ated numerically. The last entry in Table 1 is the 
difference between the logarithmically divergent 
piece which we use in the subtraction and that for 
the MS expression. 

After we get the analytic expression for the 
Feynman diagrams we listed in the above, we eval- 
uate the integrals numerically by use of VEGAS. 
Table 1 is the summary of the numerical values 
of each Feynman integral. We used a s = 0.135 
and Mb = 1.71 in lattice units. Note that Mb 
in the integral differs from that in the simulation 
(1.5) because the tadpole improvement scheme we 
employed tells us that Mb = Mo/«o- 

For the S-wave, from the relations, 

D'l at = (1 + G)D 2 + FI (15) 
I Lat = (1 + E)I, (16) 

where G = (sum of fig 2 contributions), F = (sum 
of fig 2 contributions), and E = (sum of fig 3 
contribution), we get 



(1 + G) Lat ((1 + £)(1 + G)) 
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I 



1 



Lat • 



For the P-wave, from the relations, 

A Lat ®A Lat = (l + A)A(g)A 
TLt®n at = (1 + H)T a (g) T a 
+ J A® A, 



(17) 
(18) 

(19) 
(20) 



where A = (sum of fig 5 contributions) , J = (sum 
of fig 6 contributions), and H = (sum of fig 5 
contributions), we obtain 



A (g) A 



1 



1 

fTTff) 



A Lat <g A 



Lat 



(21) 



Lat 



j 



((l + H)(l + A)) 



A Lat <g A 



Lat • 



(22) 



From [2] 

pLat 



~iLat 



we recall 

, Gi uLat 



jjLat 



0.1488(5)-^, 
1.3134(9)^, itf<* = 0.0145(6)-^, 
0.0152(3)M 6 2 #i. Thus, after putting 
all the coefficients together, we get M'^Gf 13 = 
0.147, M^Ff^ = 0.133, M^Hf s ~ = 0.01641, 
M 2 H^ S = 0.000337 in lattice units (without er- 
ror bars). With these MS matrix elements, we 
estimate the decay rates for various processes us- 
ing formulae in [1]. With a s (M b ) = 0.20, M b = 
4.7GeV, a em = 1/128, we get : For the S-wave 
case, 



T(T -> LH) 



9)(# c 2 



{1 



54# c 
-9.46(2)CV 



^-C F a 3 s {M b ) 



+4.13(17)CU - 1.161(2)^) — } 



i 

-- 36KeV 



2 V 



13a s 

47T 



C F )]2Gi 



r(T -> /+/- 



2?r 



Q 2 a 



x 16a,(M b ) 

37T 



For P-wave case, 
T( 1 P 1 -> Li?) 



^-a?(M 6 ) J ff 8 = 26KeV 



7ra*(M 6 )(C F #i + -jiff 8 ) 



T( 3 J Pi -> Li?) 



380KeV 



^a2(M 6 )ii 8 = 26KeV 



4 n f 

Ka](M b ){ — C F H 1 + -LH fi ) 

120KeV 



r( 3 i J o 77) = 67rQ 4 a2 mJ Fi = 0.26KeV 



77) 



8tt 



Q 4 a 2 mjFl = Q.070KeV 



Among these decay rates, experimental values 
exist only for electromagnetic and hadronic decay 
of the S-wave T. In this case, we can compare 
the lattice values (36KeV, 0.98KeV)) to the ex- 
perimental values (48KeV, 1.3KeV) and we find 
the agreement is reasonable. 

The remaining work is the scale setting for each 
integral listed in the above and the work is cur- 
rently in progress [4]. 
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